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Clarifying the chiral domains structure of superconducting Sr 2 Ru 04 has been a long-standing 
issue in identifying its peculiar topological superconducting state. We evaluated the critical current 
7c versus the magnetic field H of Nb/Sr 2 Ru 04 Josephson junctions, changing the junction dimension 
in expectation of that the number of domains in the junction is controlled. 7c(77) exhibits a recovery 
from inversion symmetry breaking to invariance when the dimension is reduced to several microns. 

This inversion invariant behavior indicates the disappearance of domain walls; thus, the size of a 
single domain is estimated at approximately several microns. 

PACS numbers: 74.50.-|-r, 74.70.Pq, 74.25.Sv 


Strontium ruthenate (Sr 2 Ru 04 ; SRO) [ij has long 
been studied and is now widely accepted as a spin-triplet 
superconductor. A number of experiments [2 t 6| have 
supported the pairing state of SRO as spin-triplet chi¬ 
ral p-wave with broken time-reversal symmetry. In the 
chiral p-wave symmetry, the orbital part of the pair po¬ 
tential is represented as kx In iky, which means that the 
phase of the pair potential evolves continuously by clock¬ 
wise or anticlockwise rotation in the kx-ky plane reflect¬ 
ing the finite angular momentum of the Cooper pairs. 
Thus, SRO is believed to be a typical example of a topo¬ 
logical superconductor In recent years, topologi¬ 

cal superconductivity has received considerable attention 
because Majorana fermions, which can be used for topo¬ 
logical quantum computation, are expected to emerge in 
half-quantum vortex cores or at the edges [nnii. The 
search for Majorana fermions is increasingly accelerated. 
However, the pairing symmetry of SRO is still controver¬ 
sial because some of the features peculiar to the chiral 
p-wave state, such as spontaneous magnetic fields due 
to the edge currents and chiral domains, have not been 
observed yet [Sin. 

To identify the pairing symmetry of unconventional su¬ 
perconductors, the sensitivity of the Josephson effect to 
the phase of the pair potential is quite useful. In previous 
high-temperature-superconductor experiments, the mag¬ 
netic field responses of the critical current 1^ in corner¬ 
shaped Josephson junctions and SQUIDs, whose super¬ 
conducting loop contains two interfaces with different ori¬ 
entations, have revealed the pairing symmetry to be the 
d-wave state [l7|. The same idea can essentially be ap¬ 
plied to the determination of the chiral p-wave state of 
SRO. However, prior to the detection of the chiral p- 
wave state using the corner-shaped Josephson junctions, 
the Josephson effect in SRO was not well understood; 
i.e., the conventional Fraunhofer diffraction pattern [l8l | 


of a single-boundary Josephson junction has not been 
observed yet. One of the reasons for the unconventional 
behavior is considered to be the effects of chiral domains 
and their motion during measurement. Kidwingira et al. 
reported a variety of complicated diffraction patterns in 
Pb/Cu/SRO junctions [19j. The interpretation of these 
complicated diffraction patterns is that the phase of an 
SRO crystal in a junction is spatially modulated owing 
to the existence of the chiral domains [2^ . They also re¬ 
ported several peculiar features indicating the existence 
of the chiral domain wall motion, such as a telegraph-like 
noise and a hysteresis loop in the diffraction patterns, and 
estimated the size of a single chiral domain at approxi¬ 
mately 1 pm. 

However, the size of a single domain is still a topic un¬ 
der discussion, because the estimated size is largely dis¬ 
tributed depending on the experimental probes [2ll | . The 
size larger than 50_pm was estimated by the polar Kerr 
effect experiment [2^ , while the size as small as ~400 nm 
was suggested by the scanning SQUID experiment [l6| . 
The determination of the domain size is one of the im¬ 
portant issues that could settle the pairing symmetry of 
SRO. Nelson et al. fabricated a Auln/SRO SQUID in 
which two junctions were formed at the opposite edges 
in the ab plane of an SRO crystal and reported the min¬ 
imum of the magnetic field modulation pattern in at 
zero magnetic field 0 . This result seems to suggest that 
the pairing symmetry of SRO is odd parity. On the other 
hand, Asano et al. theoretically calculated that the mod¬ 
ulation pattern can be shifted in phase by tt, depending 
on whether the number of domains in the SQUID loop is 
even or odd [i^. Since the size of the SRO crystal they 
used was on the order of millimeters which probably in¬ 
cluded a large numbers of domains, the phase shifts at 
the domain walls should be treated more explicitely. 

Here, we report the junction size dependence of the 
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FIG. 1. (a) Chiral domain structures used in the simula¬ 

tion. The two colors represent the different chiral domains 
(i.e. ±kx -|- iky) of SRO. (b) Simulation of Ic{H), taking ac¬ 
count of the self-field under domain structure A. The critical 
current Ic and the magnetic field H are normalized by Ii and 
Ho corresponding to the flux quantum <l?o, respectively. The 
dashed curve is the time-reversed version of the solid curve. 
Although the time-reversal symmetry is broken due to the 
self-field and the chiral domain, the IS is still invariant as 
long as the domain walls remain static, (c) The IS in Ic{H) 
breaks as the domain wall moves from A to B at the threshold 
field of Hthi during the field-sweep measurement. 


magnetic field H responses of the critical current Ic in 
Nb/SRO Josephson junctions. We expect that the num¬ 
ber of domains in the junction is controlled by changing 
the dimension of the junction, and correspondingly the 
diffraction pattern Ic{H) varies depending on the config¬ 
uration of domains. We focus on the inversion symmetry 
(IS) in Ic{H), which is invariant in the absence of domain 
wall motion (the details are given later). As we reduced 
the width of the junctions, the IS in Ic{H) exhibited a 
recovery from breaking to invariance at a junction width 
of several microns. This result led us to conclude that the 
size of a single domain is on the order of several microns. 


First, we illustrate with a simulation the concept be¬ 
hind our experiment. We assume a Josephson junction 
between a SRO crystal and an s-wave superconductor 
attached at a single side of SRO, and SRO to be the 
two-dimensional chiral p-wave. In the present simula¬ 
tion, we assume that the junction width is sufficiently 
smaller than the the Josephson penetration depth Aj for 
simplicity. As shown in Fig. [T](a), we employ the chiral- 
domain model in which the j/-component (parallel to the 
interface) of the pair potential keeps its phase, while the 
x-component (perpendicular to the interface) changes its 
phase by tt at the domain boundary [^. Ic{H) is eval¬ 
uated by taking the Josephson current I as the form of 
I = h cos 9 — I 2 sin 20, where Ii I 2 and 9 is the phase 
of SRO relative to that of the s-wave 23, 2J|. We calcu- 


FIG. 2. (a) Schematic of sequence used to modify the junc¬ 
tion width w. (b) SIM image of junction A at lu =10 pm. 
Nb films in yellow area were removed to fabricate the initial 
junction with w = 55 pm. After measurements at ui = 55 pm, 
the Nb films in the red areas were removed to change w to 
10 pm. (c) R-T characteristics of junction A at each w. The 
superconducting transition temperature Tc is approximately 
1.4 K at all w. (d) Typical I-V characteristic of junction A 
at to = 5.8 pm observed at T = 1.4 K [Tc ~1.41 K). 


late /c(i?)’s for positive and negative current directions 
represented by It{H) and respectively, in order 

to evaluate the symmetry of Ic{H) with respect to the 
field and the current direction. When both the chiral do¬ 
mains and the self-field are absent, Ic{H) is time-reversal 
invariant i.e., lf{H) = and Ic{H) = 

In contrast, when the effect of the self-field cannot be ne¬ 
glected, the Ic{H) calculated for chiral domain structure 
A in Fig.[T](a) exhibits broken time-reversal symmetry, as 
shown in Fig. [l](b), while the Ic{H) is still IS-invariant, 
be., I^{H) = as far as the chiral domains re¬ 

main static. On the other hand, if once the chiral domain 
wall shifts from A to B as increasing H beyond a thresh¬ 
old field ibthr at which the domain walls begin to move 
(Fig. [IKa)), the Ic{H) is modified from the black curves 
to the red curves as shown in Fig. HDc), and then Ic{H) 
is no longer invariant with respect to the IS. Therefore, 
the chiral domain motion can be detected sensitively by 
testing the IS invariance in Ic{H). We note that the 
IS is insensitive to the junction dimension, Aj, and the 
uniformity of the current. According to the previous ex¬ 
periments, chiral domain wall motion is considered to be 
excited by applied magnetic fields 0 or by a current 
flow |25l427l| . Here, we aim to test the domain wall mo¬ 
tion by reducing the junction dimension. 

Next, we move to the experimental side. We suc¬ 
ceeded in fabricating Nb/SRO Josephson junctions in 
which high supercurrent densities (as high as on the or¬ 
der of 10^ A/m^) are realized [ 2 ^. In the present work, 
the width of junctions was sequentially modified using 
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FIG. 3. Magnetic field H responses of critical current Ic{H) 
in junction A (r=0.55 K) at (a) w = 55 fim, (b) w = 
10 /rm and (c) w = 5.8 fj,m, and in junction B (r=0.4 K) 
at (d)w=10.3 fim and (e)w=5.8 pm. The applied magnetic 
field H (//c-axis) is normalized by the period for the conven¬ 
tional Fraunhofer pattern of Hq estimated for each w shown 
in the figures. I^{H) (black data) and (H) (blue data) 
are Ic{H)’s in a positive and a negative current directions, 
respectively. I^{H) (red data) is obtained from I~{H) under 
the inverse projection with respect to the current direction 
and H (represented as black arrows). The IS of the junctions 
is evaluated by the consistency between I^{H) and 
In contrast to the results at (a), (b) and (d), I~{H) is consis¬ 
tent with (H) at w = 5.8 pm in both junction A and B [(c) 
and (e)]. (f)i/-sweep-range dependence of the IS in junction 
B at ui=10.3 pm. The IS gradually recovered by reducing the 
sweep range, and Hthi is estimated at H/Hq ~ 0.25(=2.1 G). 


a focused ion beam (FIB) so that we can clarify purely 
the effect of the junction dimension. Figure UKa) shows 
the schematic illustration of the typical sequence used to 
modify the width of the junction. We fabricated junc¬ 
tions with widths w successively made narrower from w 
= 55 pm to 10 pm, and to 5.8 pm (junction A), and 
from r(;=10.3 pm to 5.8 pm (junction B). A scanning ion 
microscopy (SIM) image of junction A at re = 10 pm 
is shown in Fig. HKb). After we measured the transport 
properties of the junctions at the wider size, the junction 
width w was successively changed to the narrower size. 
The transport properties of the junctions at each w were 
measured using a standard four-terminal method down 
to approximately 0.4 K. The junctions were magnetically 
shielded by using double p-metal shields to reduce the 



FIG. 4. Hysteresis loops in Ic{H) characteristics. Magnetic 
field H was swept from zero to a negative value (black data), 
and then swept up to a positive value (blue data), and finally 
swept back to zero (red data). Diffraction patterns observed 
at u) ~ 10 pm showed hysteresis loops depending on the di¬ 
rection in which the magnetic field was swept [(a) and (c)]. 
However, the hysteresis loops disappeared in both junction A 
and B as we reduced w to 5.8 pm [(b) and (d)]. 


residual magnetic field less than 4 mG. The maximum 
critical current Ic shown later in Table |T] is almost pro¬ 
portional to w. Thus, our junctions are considered to be 
mostly uniform. The Aj’s estimated by the critical cur¬ 
rent densities are approximately 5 pm (junctin A) 
and 3 pm (junction B), respectively. Figure |2](c) shows 
the resistance-temperature (R-T) characteristics of junc¬ 
tion A at each w. Although the resistance at a normal 
state increased as w decreased, a sharp superconducting 
transition was maintained at Tc ^ 1.4 K for all w. This 
result confirms that the FIB process for modifying w did 
not damage the junction quality. As the reduction of the 
resistance at T ^2.3 K is clearly separated from the tran¬ 
sition at Tc ^1.4 K, we consider that the 3-K phase [^,0 
in the bulk SRO crystal near the junction reduced the 
resistance at T ^2.3 K. Figure |2Kd) shows a current- 
voltage (I-V) characteristic of junction A a.t w = 5.8 pm 
observed at T=1.4 K. The I-V characteristic exhibited 
a typical overdamped behavior with no hysteresis loop. 

Figure [3] shows the magnetic field H responses of the 
critical current R in junction A and B. The applied 
field (//c-axis) swept up from a negative to a positive 
value is normalized by Hq, which is the period for the 
conventional Fraunhofer pattern estimated for each w, 
Hq = $o/ [w(Asro + Ans)]) where <i>o is the flux quan¬ 
tum (20.7 G-pm^), and Asro (= 190 nm for H//c-axis) 
and ANb (= 44 nm) are the penetration depths in SRO 
and Nb, respectively [1^. The consistency or incon¬ 
sistency between I/'{H) and IV{H) determines the IS 
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invariance or breaking, respectively, where I~{H) was 
obtained from I~{H) by the inverse projection, be., 
Ic{H) = Injunction A, Ic{H) at w = 55 

tended to change irregularly, and we have not observed 
any periodic Ic{H) [Fig. [31[a)]. This behavior is reason¬ 
able because many chiral domains are considered to be 
present inside the junction area, and simultanously the 
width of w = 55/rm is much larger than Xj 5 /im. As 
the width of the junction was reduced to w = 10 /rm, the 
Ic{H) tended to exhibit a peak structure, although the 
maximum of I^{H) and the minimum of I~{H) shifted 
to a negative H as shown in Fig. jSjb). The reduction of 
the peak width from the expected value {H/Hq '^0.15) is 
attributed to the concentration of the applied magnetic 
field at the edge of the SRO crystal due to the Meissner 
effect [i^. As further reducing the junction width to w 
= 5.8 /im, the Ic{H) became rather conventinal as shown 
in Fig. EIc) ; the maximum of I^{H) and the minimum 
of Ic{H) were observed at H/Hq ^ 0. Moreover, I^iH) 
became almost consistent with indicating the re¬ 

covery of the IS invariance. This feature is quite different 
from those observed on in = 55 /xm and 10 fim junc¬ 
tions. Similar recovery of the IS has also been observed 
in junction B, be., the IS was broken at w = 10.3 /xm 
due to the difference between I^iH) and Ic{H) around 
H/Hq ~ 0 [Fig. |3l[d)]. However, a conventional pattern 
with the IS invariance was recovered at w = 5.8 /xm, as 
shown in Fig.[3l[e). In order to confirm the validity of the 
chiral-domain model, we further evaluate the threshold 
field ibthr in junction B. Figure [3l[f) shows the magnetic 
field-sweep-range dependence of Ic{H). The data of the 
largest sweep range in Fig. [HJf) is same as that shown 
in Fig. Eld). The gradual recovery of the IS invariance 
by reducing the field-sweep range means the suppression 
of chiral domain motion in smaller field. The complete 
recovery of the IS in the lowest curve indicates that the 
Hthr/Ho ^ 0.25 (ibthr=2.1 G) in junction B. Therefore, 
the recovery of the IS invariance at w = 5.8 /xm under 
the field-sweep range of ±14 G, which is far larger than 
ibthr, shown in Fig. He) suggests the absence of domain 
walls inside the junction area. 


In addition to the IS, we detected the disappearance 
of a hysteresis loop in both junction A and B as xc was 
reduced. The magnetic field was swept from zero to a 
negative value, and then swept up to a positive value, and 
finally swept back to zero. At xc ^ 10 /xm, we observed a 
hysteresis loop depending on the direction in which the 
magnetic field was swept [Fig.HDa) and (c)]. In Fig.HDa), 
the maximum of (H) and the minimum of I~ (H) shift 
in the direction in which the magnetic field was swept. 
In Fig. |4Dc), the magnitude of the maximum I^{H) and 
the minimum I~(H) changes depending on the sweep 
direction. These hysteresis loops disappeared at xc = 
5.8 /xm in the both junction A and B, as shown in Fig. 
SDb) and (d). Accepting that the origin of the hysteresis 
loop is the chiral domain wall motion [3113 , the chiral 


TABLE I. Summary of junction size dependence of magnetic 
field H responses of critical current 7c. The value of maximum 
/c is averaged over several cooling cycles, and A/c/7c, where 
A/c is the standard deviation, is estimated. 


Junction 

w [/xm] 

Ic [mA] 

AIc/Ic [%| 

IdH) 

IS 

A 

(T=0.55 K) 

55 

9.33 

8.39 

Random 

X 

10 

1.55 

1.09 

Hysteresis 

X 

5.8 

0.94 

0.58 

Conventional 

o 

B 

(r=0.4 K) 

10.3 

11.71 

2.64 

Hysteresis 

X 

5.8 

7.96 

1.33 

Conventional 

o 


domains are considered to be present at xc ^ 10 /xm, while 
they disappear at the junctions of 5.8 /xm. 

Table U summarizes the results of the junction size 
dependence of Ic- In addition to the above mentioned 
features, we discuss the distribution of Ic estimated by 
IXlc/Ici where A7c is the standard deviation of Ic over 
several cooling cycles. Assuming that chiral domain tex¬ 
tures are expected to be inequivalent in each cooling cy¬ 
cle, the variation in Ic reflects the presence of multiple 
numbers of chiral domains and the variation of their con¬ 
figuration. Thus, the tendency that AJc/7c decreases as 
w is reduced reflects that the number of chiral domains 
decreases as w is reduced. Putting together the IS invari¬ 
ance and the lack of the hysteresis loop in both junction A 
and B as xn is reduced, we conclude that the detected size 
dependence of Ic{II) is governed by the chiral domains 
and their motion, and that the size of a single chiral do¬ 
main is estimated on the order of several microns. 

The domain size is almost consistent with those es¬ 
timated by several results using the 3-K phase 
while it is somewhat larger than ~I /xm estimated by 
Kidwingira et al. [13 ■ We speculate that the relatively 
larger domain can be induced by s-wave Nb films whose 
superconducting transition temperature Tc (~6.5 K) is 
higher than that of SRO. In our Josephson junctions, 
the contact between the Nb films and SRO realizes a 
high current density [23 compared to that of other junc¬ 
tions using the I.5-K phase. Thus, the phase of SRO 
was locked to that of the Nb film through the Joseph¬ 
son coupling, which probably results in the creation of 
the relatively large size of domains. The effect of the 
phase lock can be checked by using s-wave superconduc¬ 
tors whose Tc is lower than that of SRO, such as Al (typi¬ 
cal Tc ~1.2 K). This phase-lock effect also might be asso¬ 
ciated with the difference in IdH) over multiple cooling 
cycles. The IciH) of junction B {w=5.8 /xm) was com¬ 
pletely stable against the several cooling cycles, whereas 
that of junction A (x(;=5.8 /xm) showed a different pattern 
for the rare occasion. Since the current density of junc¬ 
tion B is about one order of magnitude higher than that 
of junction A, it is reasonable to conclude that the stable 
7c(77) in junction B is also attributed to the phase-lock 
effect. 
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Although the chiral p-wave symmetry of SRO has been 
assumed throughout this paper, the present result does 
not exclude the possibility of the helical p-wave symme¬ 
try that is another candidate of the paring symmetry of 
SRO We believe that the analyses are mostly un¬ 
changed if we consider the presence of helical domains 
instead of chiral domains. For the analysis based on the 
helical domain models, theoretical calculation performed 
on the helical domain boundaries is required. 

In summary, we measured the junction size dependence 
of the magnetic field H responses of the critical current 
Ic in Nb/Sr 2 Ru 04 Josephson junctions and tested the in¬ 
version symmetry (IS) invariance of Ic{H). The IS exhib¬ 
ited a recovery from breaking to invariance at a junction 
width of several microns. This inversion invariant recov¬ 
ery indicates the absence of chiral domain wall motion 
and led us to conclude that the size of a single chiral do¬ 
main is on the order of several microns. These results will 
open the possibility that the internal phase of Sr 2 Ru 04 
can be identified by using corner-shaped junctions with 
a size smaller than several microns in the future. 

We are grateful to M. Koyanagi for the fabrication of 
the junctions. We thank S. Yonezawa, and M. S. An¬ 
war for fruitful discussions. This work was supported by 
MEXT KAKENHI Grant Numbers 22103002,15H05852 
and 15H05853, and by Grants-in-Aid for Scientific Re¬ 
search (No. 60443181 and No. 15H03686) from the Japan 
Society for the Promotion of Science, Japan. 


* current affiliation: Department of Applied Physics, Os¬ 
aka University, Suita 565-0871, Japan 

[1] Y. Maeno, H. Hashimoto, K. Yoshida, S. Nishizaki, 
T. Fujita, J. G. Bednorz, and F. Lichtenberg, Nature 
372, 532 (1994). 

[2] K. Ishida, H. Mukuda, Y. Kitaoka, K. Asayama, Z. Q. 
Mao, Y. Mori, and Y. Maeno, Nature 396, 658 (1998). 

[3] G. M. Luke, Y. Fudamoto, K. M. Kojima, M. I. Larkin, 

J. Merrin, B. Nachumi, Y. J. Uemura, Y. Maeno, Z. Q. 
Mao, Y. Mori, H. Nakamura, and M. Sigrist, Nature 
394, 558 (1998). 

[4] K. D. Nelson, Z. Q. Mao, Y. Maeno, and Y. Liu, Science 
306, 1151 (2004). 

[5] Y. Maeno, S. Kittaka, T. Nomura, S. Yonezawa, and 

K. Ishida, J. Phys. Soc. Jpn. 81, 011009 (2012). 

[6] A. P. Mackenzie and Y. Maeno, 
Rev. Mod. Phys. 75, 657 (2003) 

[7] S. Kashiwaya, H. Kashiwaya, H. Kambara, T. Fu- 


ruta, H. Yaguchi, Y. Tanaka, and Y. Maeno, 
Phys. Rev. Lett. 107, 077003 (2011) 

[8] T. Nakamura, R. Nakagawa, T. Yamagishi, 
T. Terashima, S. Yonezawa, M. Sigrist, and Y. Maeno, 
Phys. Rev. B 84, 060512 (2011) 

[9] T. Nakamura, T. Sumi, S. Yonezawa, T. Terashima, 

M. Sigrist, H. Kaneyasu, and Y. Maeno, J. Phys. Soc. 
Jpn. 81, 064708 (2012). 

[10] J. Jang, D. G. Ferguson, V. Vakaryuk, R. Budakian, S. B. 
Chung, P. M. Goldbart, and Y. Maeno, Science 331, 186 
( 2011 ). 

[11] N. Read and D. Green, Phys. Rev. B 61, 10267 (2000) 

[12] X.-L. Qi and S.-C. Zhang, 

Rev. Mod. Phys. 83, 1057 (2011) 

[13] J. Alicea, Rep. Prog. Phys. 75, 076501 (2012). 

[14] Y. Tanaka, M. Sato, and N. Nagaosa, J. Phys. Soc. Jpn. 
81, 011013 (2012). 

[15] J. R. Kirtley, C. Kallin, C. W. Hicks, E.-A. Kim, 
Y. Liu, K. A. Moler, Y. Maeno, and K. D. Nelson, 
Phys. Rev. B 76, 014526 (2007) 

[16] C. W. Hicks, J. R. Kirtley, T. M. Lippman, 

N. C. Koshnick, M. E. Huber, Y. Maeno, 
W. M. Yuhasz, M. B. Maple, and K. A. Moler, 
Phys. Rev. B 81, 214501 (2010) 

[17] D. J. Van Harlingen, Rev. Mod. Phys. 67, 515 (1995) 

[18] A. Barone and G. Paterno, Physics and Application of 
the Josephson Effect (Wiley, New York, 1982) Chap. 4 - 
5. 

[19] F. Kidwingira, J. D. Strand, D. J. Harlingen, and 
Y. Maeno, Science 314, 1267 (2006). 

[20] A. Bouhon and M. Sigrist, New J. Phys. 12, 043031 

( 2010 ). 

[21] C. Kallin and A. J. Berlinsky, J. Phys.: Condens. Matter 
21, 164210 (2009). 

[22] J. Xia, Y. Maeno, P. T. Beyersdorf, M. M. Fejer, and 
A. Kapitulnik, Phys. Rev. Lett. 97, 167002 (2006) 

[23] Y. Asano, Y. Tanaka, M. Sigrist, and S. Kashiwaya, 
Phys. Rev. B 71, 214501 (2005) 

[24] Y. Asano, Y. Tanaka, M. Sigrist, and S. Kashiwaya, 
Phys. Rev. B 67, 184505 (2003) 

[25] H. Kambara, S. Kashiwaya, H. Yaguchi, 

Y. Asano, Y. Tanaka, and Y. Maeno, 

Phys. Rev. Lett. 101, 267003 (2008) 

[26] H. Kambara, T. Matsumoto, H. Kashiwaya, S. Kashi¬ 
waya, H. Yaguchi, Y. Asano, Y. Tanaka, and Y. Maeno, 
J. Phys. Soc. Jpn. 79, 074708 (2010). 

[27] M. S. Anwar, T. Nakamura, S. Yonezawa, M. Yakabe, 
R. Ishiguro, H. Takayanagi, and Y. Maeno, Sci. Rep. 3, 
2480 (2013). 

[28] K. Saitoh, S. Kashiwaya, H. Kashiwaya, M. Koyanagi, 
Y. Mawatari, Y. Tanaka, and Y. Maeno, Appl. Phys. 
Express 5, 113101 (2012). 





